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(N Abstract 

^P 1 A scattering problem (or more precisely, a transmission-reflection problem) of linearized excitations in the presence of a dark 
soliton is considered in a one-dimensional nonlinear Schrodinger system with a general nonlinearity: \d t <p = -d 2 cf> + F(\<p\ 2 )cf>. 



If 



the system is interpreted as a Bose-Einstein condensate, the linearized excitation is a Bogoliubov phonon, and the linearized 
C*~) equation is the Bogoliubov equation. We exactly prove that the perfect transmission of the zero-energy phonon is suppressed at a 
7 ~ ' critical state determined by Barashenkov's stability criterion [Phys. Rev. Lett. 77, (1996) 1193.], and near the critical state, the 

| ! energy-dependence of the reflection coefficient shows a saddle-node type scaling law. The analytical results are well supported by 

numerical calculation for cubic-quintic nonlinearity. Our result gives an exact example of scaling laws of saddle-node bifurcation in 



time-reversible Hamiltonian systems. As a by-product of the proof, we also give all exact zero-energy solutions of the Bogoliubov 



I equation and their finite energy extension. 
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& 1. Introduction 
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In this paper, we solve a scattering problem of lin- 
earized excitations in the presence of a dark soliton in one- 
dimensional(lD) nonlinear Schrodinger (NLS) equation with a 



CJ general nonlinearity: 
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(1.1) 



and discuss the physical and mathematical significance of our 
results. For a schematic picture, see Fig. Q] The precise math- 
ematical definition of the problem will be given in Sec. [2] If 
we regard the system as a Bose-Einstein condensate(BEC), the 
linearized excitation is a Bogoliubov phonon, so the problem 
can be also called a soliton-phonon scattering problem, as this 
paper entitled. 

The NLS equation dl.ll ) has a great number of applications 
in nonlinear optics, superconductors, magnetism, BECs, and so 
on. Particularly, much attention has been focused on the ex- 
perimental realizations of BECs in laser-cooled ultracold atoms 
for more than a decade, because of high-controllability of the 
system parameters. By using elongated laser beams, low- 
dimensional systems are realized, and a dark soliton can be 
created via the phase imprinting methodJU. The Bogoliubov 
theory is also well confirmed QH, 01- 

It is known that ID NLS with a cubic nonlinearity is com- 
pletely integrablej5l0]. Because of integrability, the linearized 
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equation is also solved exactly 10, 01 . an d tne phonon excita- 
tions are shown to be completely reflectionless against a soli- 
ton for any excitation energy. Thus, the problem is trivial in 
this case. However, when the nonlinear term is generalized, 
the phonon has a finite reflection coefficient in general. It is 
worthy to note that the soliton decay dynamics in the laser- 
trapped quasi- ID BEC has been well explained by the quintic 
term, which a ppe ars as a second-order perturbation of the trap- 
ping effect! 9, llOUllll . and yields the frictional force between 
thermal excitation clouds and soli tons Thus, knowing 

the scattering properties between solitons and linearized exci- 
tations is essential to understand and control the transport of 
solitons, that is, the transport of stable wave packets. We also 
mention that the theory of nonpolynomial NL S eq uation is for- 
mulated to describe the confinement effect B12lll3Lll4ll5ll . The 
quintic NLS also appears in an effective mean-field description 
of the Tonks-Girardeau gas lfl^l . 

The NLS equation with an integrability breaking factor is 
also interesting from the viewpoint of an infinite-dimensional 
dynamical system and the bifurcation theory. When the poten- 
tial barrier is added in the cubic NLS equation, there exist stable 
and unstable stationary supercurrent-flowing solutions 11 71 1 1 811 . 
if the condensate velocity does not exceed a certain critical 
value. Near the critical point, which separates the stable branch 
and the unstable branch, it is known that many physical quanti- 
ties obey saddle-no de ty pe scaling, such as an emission period 
of dark solitons lflcl Il9tl. an eigenvalue of a growing mode for 
the unstable solution] 19], and a transmission coefficient of lin- 
earized excitations! 8, 20, 21]. It is quite nontrivial that the time- 
reversible Hamiltonian system exhibits the scaling behaviors of 
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saddle-node type, since this bifurcation is normally understood 
to emerge in time-irreversible phenomena. However, it is not 
easy to prove these properties analytically or exactly, because 
of the infinite dimensional character of the system. 

On the other hand, as another way to break the integra- 
bility, one can consider the generalization of the nonlinear- 
ity, that is what we will consider in the present paper. When 
the nonlinear term includes a competing interaction, the dark 
soliton is no longer always stable. One typical example of 
an unstable dark soliton is a "bubble" in a cubic-quintic NLS 
(CQNLS) systemJUIH. (See also Q.) The most general 
criterion for the stability of the dark soliton has been shown by 
B arashenkov 12511 . and the critical velocity of the soliton is de- 
termined by dP/dv = 0, where v is a velocity of the soliton 
and P is a renormalized momentum. The existence of the crit- 
ical velocity lower than the sound velocity and the separation 
of stable and unstable regions are similar to the phenomena of 
superflows against a potential barrier. Therefore, we can expect 
some scaling behavior near the critical state. Furthermore, in 
the present case, the preserved translational symmetry of the 
fundamental equation makes it possible to access the problem 
analytically. 

In this paper, we solve the scattering problem of linearized 
phonon excitations, and exactly show the following: (i) At 
the critical state determined by Barashenkov's criterion [25], 
the transparency of the zero-energy phonon is suppressed, and 
only partial transmission occurs, (ii) Near the critical state, the 
energy-dependence of the reflection coefficient of low-energy 
phonons shows saddle-node scaling behavior, regarding the 
renormalized momentum as a parameter of a normal form of 
saddle-node bifurcation. The obtained analytical results are 
well confirmed by comparison with the numerical results of 
CQNLS equation. Our result gives an exact example of scaling 
laws of saddle-node bifurcation in time-reversible Hamiltonian 
systems. The proof is based on the exact low-energy expan- 
sion of the solution of the linearized equation. Since the exact 
zero-energy solutions given in this paper are quite general, we 
believe that our method will also be useful to derive other low- 
energy physical properties. 

The organization of the paper is as follows. In Sec. [2j we in- 
troduce fundamental equations and see the fundamental prop- 
erties. The definition of the transmission-reflection problem is 
also given. In Sec. [3] we give a main result and verify it by nu- 
merical study of CQNLS equation. Sections|4]and[5]are devoted 
to the proof of main results. Discussions, future perspectives, 
and conclusions are given in Sec. [6] Some mathematically tech- 
nical formulae are treated in Appendices. 



2. Fundamental Equations and Definition of the Problem 

2.1. Fundamental equations 

We begin with the NLS equation with a general nonlinearity 

i3,0 = -c^ + F(|0| 2 )0. (2.1) 




Figure 1: A schematic picture of the problem that we consider in this paper. 
p represents a half of the velocity of the dark soliton. The problem is always 
considered in the comoving frame of the soliton. e'* i r is an incident wave of a 
linearized excitation, fe'*' r is a transmitted wave, and re'* 2 > is a reflected wave. 
For more detailed definitions of each quantity, see Sec. [2] 



Here, F(p) is a real-valued function such that F(0) = 0. The 
energy functional (Hamiltonian) which yields this equation is 

H= Jdx(|5,0| 2 + C/(|0| 2 )), (2.2) 



where 



U(p) 



P ' dp'F(p'). 
Jo 



(2.3) 



Letting (f> = (p + 6(p in Eq. ( 12. It . and discarding higher order 
terms of 6(p, one obtains the following linearized equation: 



id,6cf, = [-dl + F(\<f>\ 2 ) + \<t>\ 2 F'(\<f>\ 2 )] 6<f> + cf> 2 F'(\cl>\ 2 )6</>*. 



Writing 5(f) - u and -5<p* = v, one obtains 



id, 



(2.4) 



(2.5) 



where X. is a 2 x 2 matrix operator whose components are 

Ai = --C22 = ~dl + F(W 2 ) + I0I 2 ^'(I0I 2 ), (2.6) 
A2 = --C21 = -<P 2 F'(\<f>\ 2 ). (2.7) 

We use the notation (u, v) since it is commonly used by con- 
densed matter physicists. If we interpret the system as BEC, 
this equation is the Bogoliubov equation which describes the 
Bogoliubov phonon (or Bogoliubov quasiparticle) IT261] . (As a 
review or a textbook, see, e.g., [27, 28, 29].) For this reason, 
henceforth, we call the condensate wavefunction or the order 
parameter, and (m, v) the Bogoliubov (quasiparticle) wavefunc- 
tion, though the NLS equation itself has more applications in 
various fields. 

Henceforth we mainly consider the stationary (i.e., time- 
independent) problem. The stationary NLS equation with 
chemical potential [i is obtained by setting 4>(x, t) = <p(x)e~ 1,ir : 



(- i u-^ + J F(|0| 2 ))0 = O. 



(2.8) 



As will be seen, the value of fi is fixed by the asymptotic form 
of (p . The stationary Bogoliubov equation with eigenenergy 



2 



e is obtained by setting u(x, t) = u(x)e and v(x, t) - 



with 



-n 

p. 



(2.9) 



(2.10) 



2.2. Bogoliubov phonons in a uniform condensate 

Let us derive the dispersion relation (the energy-momentum 
relation) of Bogoliubov phonons when the condensate is flow- 
ing uniformly: (p(x) = -\fp^e l(j,x+lf \ In order for this <p{x) to be 
the solution of Eq. I I2.8I ). the chemical potential must be 

fX = p 2 + F(p x ). (2.11) 

The four solutions of Bogoliubov equation ( 12.91 1 are given by 



"fort := v r 



(2.12) 



where i — 1, 2, 3, and 4, and the wavenumber k,s are the roots 
of the following quartic equation: 

(e - 2kp) 2 = k 2 ^ 2 + 2p x F'( Poo )). (2.13) 



Equation ( 12.131 ) gives the dispersion relation, and from this dis- 
persion one can see that a half of the Landau's critical velocity 
(or a half of the sound wave velocity) is given by 



PL = 



PooF'ipoc) 



(2.14) 



(Note that the sound wave velocity is not /?l but 2p^\ see the 
next subsection.) Since p^ must be real, in order for the uniform 
condensate to be stable, F'(p<x,) > must hold. The coefficients 
w; and v, can be, e.g., chosen as follows: 



_ e-2pL + k] 

Ui= f + 2p 2 L ' 

_ f e-2pk.-k] 

Vi= f 



(2.15) 



(2.16) 



When e > and — Pl < p < p^, the quartic equation ( 12. 1 3b 
has one real positive root, one real negative root, and two com- 
plex roots conjugate to each other. We call a real positive (neg- 
ative) root k\ (kx), and a complex root with positive (negative) 
imaginary part kj, (£4). The low-energy expansions of them are 
given by 

e 



h = 



+ 0(6% 



2(p + Pl) 

k 2 = + Oie 3 ), 

2(P ~ Pl) 



h = 2i yjpl - p 2 + 



pe 



2(p 2 L - p 2 ) 



0(e 2 ), 



p z + 



2(Pl ~ P 2 ) 



(2.17) 
(2.18) 
(2.19) 
(2.20) 



w\ and W2 are plane wave solutions propagating in the positive 
and negative directions, respectively. W3 and w\ are exponen- 
tially divergent unphysical solutions. 

2.3. Dark soliton solution in comoving frame 

Let us consider the dark soliton solution of stationary NLS 
Eq. ( 12.81 1 in the comoving frame of the soliton. In this coor- 
dinate, the soliton is static but the surrounding condensate is 
flowing. Let us seek the solution with the asymptotic form 



<p(x — > ±00) = yfpZe 



iQ«±§) 



(2.21) 



It should be noted that the velocity of the soliton is not —p but 
-2p, because the Galilean covariance of NLS equation leads to 
the following property: 



<p(x, t) is a solution. 
<-> 0(jc, t, a) = cf>(x + 2at, t)eT lax eT la ~' is a solution 



(2.22) 



So, if one has the solution in the form <p(x, f) = e~ lf "e lpx f(x), the 
corresponding soliton-moving solution is given by 0(jc, t, p) = 
e -i(f-p-)i y( x + 2pt). However, for brevity, we sometimes call p 
"velocity", ignoring the difference of twice factor. 

From the conservation laws of mass and momentum, one can 
immediately find two integration constants: 



J = 



0*<p x - (j)<t>* x 
2i 



j m = \4> x \ 2 + p\<P\ 2 - U(\cf>\ 2 ). 



(2.23) 
(2.24) 



Here j is a mass current density and j,„ is a momentum current 
density. Let us write the density and the phase of the condensate 
as <p = ~\/pe lS . Taking account of the asymptotic form ( I2.21I I. 
the chemical potential p becomes the same as 1 12.1 II ). and the 
above constants are determined as 



J = PooP, 

j,„ = 2p oa p 2 + pco-F(pco) - UQjoc). 

The conservation laws are then rewritten as 

Pood* 



j P oop r x , 

s x = - = — <-> s = p \ 

P P Jo 



(2.25) 
(2.26) 



(2.27) 



(Px) 2 



= -p 2 {Poo -pf +p[U(p)~ U(poc) - (P-Pcc)U'(Pco)] ■ 

(2.28) 



Thus, one can at least obtain the formal solution 

dp 



+2(x - xq) 



I Jrm 



y/R.H.S. ofEq. QJgb 



(2.29) 



even though it is not easy in general to carry out this integration 
and obtain the solution in closed form "p(x) = . . . ". Hence- 
forth we do not need this formal solution, but we assume the 
existence of a dark soliton solution which has no singularity 
and satisfies the asymptotic condition ( 12.211 ). 
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From Eq. (12.27b . The phase shift 5 in Eq. ( 12.21b can be writ- 
ten down explicitly: 



P 



(2.30) 



We also introduce the symbol for the particle number of the 
dark soliton for later convenience: 



•J —• 



N := &x(p- p x ) <0. 



(2.31) 



2.4. Barashenkov's criterion 



The stability of the soliton is described by the following 
renormalized momentum|25]: 



2i 



(2.32) 



Here (j>(x, t) is the dark soliton solution in the frame where the 
surrounding condensate is at rest and the soliton is moving. 
Writing the soliton velocity v, the stability criterion for the dark 
solitons is expressed by dP/dv < 0. 

We can rewrite the above integral by the density profile p(x): 



T 

= -pN - poo 6. 



P= - P \ civ:/>, 



(2.33) 



Here remember that the soliton velocity is given by v = —2p, 
as stated in the preceding subsection. The stability condition is 
rewritten as d(—P)/dp < 0. 

2.5. Definition of the scattering problem 

In this subsection, we define the transmission and reflection 
problem of Bogoliubov phonons shown in Fig. [TJ Since the lin- 
earized equation does not satisfy simple particle number con- 
servation, we must define transmission and reflection coeffi- 
cients via the conservation of excitation energy. The conser- 
vation of excitation energy corresponds to the constancy of the 
following Wronskian: 1201 12 ill : 



W : 



*d x u — ud x u* + v*d x v — vd x v* . 



(2.34) 



Let us assume that the asymptotic form of the condensate wave- 
function <p(x) is given by Eq. ( 12.21b . In this situation, suffi- 
ciently far from the origin, the Bogoliubov equations have the 
plane wave (and exponentially decaying/diverging) solutions 
given by Eq. ( 12.12b with <p = ±1, 

The solution of the scattering problem is defined by the one 
that has the following asymptotic formJi, 2(| 21]: 



wi (*,-§) + rw 2 (*,-§) fx->-oo) 
tWj (x, +§) (x -> +00). 



(2.35) 



Here the exponentially decaying waves, which are Wzt(x, — f) 
-00 and W3(x, +§ ) in x — > +00, can be also included. 



in x 



However, they are irrelevant in the definition of transmission 
and reflection coefficients. The calculation of W shows that 



=|f| 2 [(ki + pMif + (*i - p)|v!| 2 ], (2.36) 

(2.37) 



2i 

W(-oo) 



2i 



^ 1 +p)|S 1 | 2 + (^ 1 - j p)|v 1 | 2 



\r\ 2 [(k2+pm 2 \ 2 + (k2-p)\v 2 \ 2 ]. 



Since W(+oo) = W(-oo), the transmission coefficient T and the 
reflection coefficient R are naturally defined as 



T = \t\\ 



(2.38) 



(ki + p)\ui\ 2 + (ki -p)\vi\ 2 

By this definition, T + R = 1 always holds. If one chooses the 
normalization S; and v,- as Eq. ( 12.15b and ( 12.16b . one can show 

(-k 2 - p)\u 2 \ 2 + {-h + p)\v 2 \ 2 = j _ pp L £ 2 + _ 4- 



(*i + P )\u i \ 2 + (k l - P )\n\ 2 



2( P 2 -p 2 L ) 3 



(2.40) 



So, if one is only interested in the leading order, one can ap- 
proximate R a \r\ 2 . 



3. Summary of Main Result and Numerical Verification 

In this section we present the main results of this paper and 
verify them by numerical study of CQNLS equation. The proof 
will be given in Secs.|4]and|5] 

3.1. Main result 

In the scattering problem of linearized excitations defined in 
Subsec. 12.51 the amplitude of the reflected component r in Eq. 
( 12.35b is given by the following Pade approximant-like form: 



-ie(d + diP p ) + 0(e 2 ) 
aP p -\e{b + biP p ) + 0(e 2 )' 



Here X p :- dX/dp and 



a = 4p L p oa , 

b = (N + pN p ) 2 + ( PL N p ) 2 , 
d = (N + pN p ) 2 - (p L N p ) 2 , 
P 2 + P 2 - 



pI-p 2 



with 



d\ - N + N 



dN dN 
N := p— -poo- — ■ 

op dpco 



(3.1) 

(3.2) 
(3.3) 
(3.4) 

(3.5) 

(3.6) 

(3.7) 
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From ( 13.11 ). the energy dependence of the reflection coefficient 
R 1 12.391 ) becomes 

i2 



R 



Id + diPpY , , 



(3.8) 



(P P = 0). 



When P p ^ 0, the zero-energy phonon transmits perfectly: 
lim e ^o R — 0. On the other hand, when the soliton velocity 
reaches a critical value, i.e. P p = 0, this perfect transmis- 
sion disappears. We note that the rational form expression d3.ll ) 
makes it possible to unify the description of low-energy behav- 
iors in both critical and non-critical cases. If we use a simple 
Taylor series, the singular behavior at the critical velocity state 
cannot be expressed. 

The above (13.8b is one good result valid for any soliton ve- 
locity, even if it is far from the critical state. However, when the 
velocity comes close to the critical value, i.e., P p comes close to 
zero, we can derive a more powerful scaling law as below. Let 
us assume that coefficients of e" (n > 2) in d3.lt are all finite in 
the limit P p — > 0. and take the limit e — ♦ and P p — > with a 
constraint e/P p = fix. We then obtain 



-id(e/P p ) 



(3.9) 



a-ib(e/P p )' 

and in the same limit, the universal form of reflection coefficient 

d\ejP p f 



lim R - , 

e-KJ^-K), a 2 + b 2 (e/P p ) 2 



(3.10) 



follows. Here, the values of a, b, and d in the critical state must 
be substituted when we use Eq. (13. 10b . We remark that Eq. 
d3.10l ) contains only ^-derivatives, whereas the expression be- 
fore taking the scaling limit contains p,x,-derivatives in addition 
to /^-derivatives. 

Let p c be a critical velocity of the dark soliton, i.e., P'(p c ) = 
0. The expansion of P near p = p c gives 



1/2 



P(P) - P(Pc) + 2 P "(Pc)(P - Pcf +"', 
-» 1^1 = \P'(P)\ * \2P"(p c )(P(p) - P( Pc ))\ 

Therefore we obtain 

e e 

]Ppi * \2P"(p c )(P(p)-P(pM /2 ' 
This is an expected scaling behavior from the normal form of 
saddle-node bifurcation! 3(J HI], if we regard the renormalized 



(3.11) 
(3.12) 

(3.13) 



momentum P as a parameter of normal form. 

3.2. Comparison with Numerical Results in CQNLS System 

In this subsection, we numerically verify the analytical re- 
sults of the preceding subsection in the CQNLS system. We 
first derive the expressions for the dark soliton solution and the 
renormalized momentum in Subsec. 13.2.11 and solve the scat- 
tering problems of linearized excitations for (i) the purely cubic 
case in Subsec. 13.2.21 (ii) the purely quintic case in Subsec. 
13.2.31 and (iii) the case where a non-trivial critical velocity ex- 
ists in Subsec. 13.2.41 



3.2.1. Dark soliton solution and renormalized momentum 
In the CQNLS system, the nonlinear term is defined by 

U(p) = flip 2 + a 2 p\ (3.14) 
F(p) = U'{p) = 2a iP + 3fl 2 p 2 , (3.15) 
and the NLS equation (12. 1 ft has the cubic-quintic nonlinearity: 
id,(p = -d\(f> + 2fl!|0|V + 3fl 2 |0lV- (3-16) 

The stationary linearized equation, i.e., the stationary Bogoli- 
ubov equation d2.9b is given by 



-d 2 - H + 4ai|^| 2 + 9a 2 |0| 4 -2ai<£ 1 - 6a 2 \cf>\ 



2a\(, 

u\ 

v 



+ 6a 2 \cf>\ z 



dl+P 



f - 6a 2 \cf>\ 2 4> 2 
4 ai \4>\ 2 -9a 2 \4>\ 4 



(3.17) 



It is known that a bubble and unstable dark solitons appear 
when «i < and a 2 > 1221 12311 . This case is considered in 
Subsec. [3.2.41 As shown in 1 23] , when the soliton velocity is 
smaller than the critical value, a small perturbation induces "nu- 
cleation dynamics", and the soliton cannot preserve its shape 
any more. So, this instability is not convective but absolute. 
The Landau velocity d2. 14b is given by 



PL = VPc°( fl l + 3fl2Poo). 



(3.18) 



The necessary condition a\ + 3a 2 pca > follows for a uniform 
condensate to be stable. The dark soliton solution is given by 



(x,p,poo) = e 



kpq + ipipoo - po) tanfiKx 
-^Po(k 2 - fl2(Poo - po) 2 tanh 2 kx) 



(3.19) 



with 



K = ^P 2 L~P 2 > 

po = p(x = 0) = 



(3.20) 

-(2a 2 p m + a\) + V(2a 2 Poo + fli) 2 + 4a 2 p 2 



2a 2 



(3.21) 



See Appendix G for a detailed derivation. Since k and po are 
the functions of (p,p<x,), the dark soliton solution has two pa- 
rameters (p,pco). From (13.21b . 



limpo 

p->0 







(la-ipoo +a\ > 0) 



2~ \2a 2 pao + ail (2a 2 pao + a x < 0), 



(3.22) 



so the bubble (= a non-topological dark soliton) appears when 
Pco < -a\l(2a 2 ). A particle number of soliton N and a phase 
difference 6 are calculated as 



N : 



V«2 



6-2 tan 



tanh 



pifioc - po) 



V"2(P°° - Po) 



Po* 



(3.23) 
(3.24) 



From them we can calculate the renormalized momentum —P = 
pN + pood. An example is shown in Fig. [2] The case where the 
unstable region exists is analyzed in Subsec. I3.2.4l in detail. 
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Figure 2: (Color online) Plot of renormalized momentum -P/p„ =5 + pN/p m 
in CQNLS. Here we set (a\,ai) = (-1, 1). The values of p^ of each curve are 
set, from top to bottom, p„ = 0.55, 0.52, 0.502, 0.5, 0.498, 0.48, and 0.45, 
respectively. The dark solitons are stable in the regions of the solid lines, while 
unstable in the regions of the dashed lines. The critical points are marked by 
black dots. The unstable soliton appears when p x < -axlilai) = 0.5. 
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Figure 3: (Color online) Energy-dependence of reflection coefficient R of lin- 
earized excitations for various soliton velocities in the purely quintic system. 
Here we set (a\,ai) = (0, 1) and p^ = 1. The Landau velocity is given by 
PL = V3poo = V3. Parabolic curves represent the theoretical approximate 
expression )3.8t with )3.34t . 



3.2.2. Purely cubic case 

As a first example, let us consider the case a 2 = 0, i.e., the 
nonlinearity is purely cubic. As mentioned in the Introduction, 
the NLS equation is integrable in this case and the Bogoliubov 
phonons are reflectionless for any energy. Let us see that our 
analytical result Eq. ( 13.81 1 is consistent with these known facts. 

Without loss of generality, we can set a\ = 1. The dark 
soliton solution ( 13. 191 ) is then reduced to 



<p — e ipx (p + ik tanh kx) , k = -^Ooo - p 2 - (3.25) 

The exact solution of the linearized equation, i.e., the Bogoli- 
ubov equation is given by the squared Jost solutionj7l[M]: 

M = e %+^Ltanh/cx+^ + — ) , (3.26) 
{ 2 2kjj 



v = e l(kj p)jc |ktanh^ + 



2 2kj 



(3.27) 



where kjs are given by the roots of the dispersion relation Eq. 
( T2~T3l > with F'(p m ) = 2. 

From the above explicit expression, it is obvious that the 
phonons are reflectionless. Therefore, the coefficient of e 2 in 
Eq. (13.81 ) must vanish. Let us check it. For the cubic case, it 
follows that 



N = -2k, 6 = 2 tan" 



(3.28) 



by taking the limit a 2 -> of Eqs. ( EHBl and l U54b . With the 
use of them, we can show 

P p = -N - pN p - Poo 6 p = 4k, (3.29) 
d = 4(k 2 - 3p 2 ), di = ~ K +3p . (3.30) 

K 

Thus we obtain d + d\P p = 0, as expected. We also note that 
the soliton is always stable since —P p < for all velocities. 



It is also possible to discuss the reflection properties when 
the quintic term is small by expanding Eqs ( 13.231 ) and 0.241 ) 
with respect to a 2 , but the expression is not so simple. In this 
case, one can derive an approximate formula valid not only for 
small energy but for arbitrary energy by the method given in 
Refs. Iffi. 

3.2.3. Purely quintic case 

Next, we treat the purely (self-defocusing) quintic case. As 
already mentioned, the quintic NLS equation is known to de- 
scribe the dynamics of the Tonks-Girardeau gas lfl6ll . 

Without loss of generality, we can set a\ = and a 2 — 1. 
Though we can also set = 1 , we keep it for a moment be- 
cause we need the differentiation of p,x, to calculate the reflec- 
tion coefficient ( 13.81 ). Eqs. (13.23b and (13.24b are reduced to 



N : 



• tanh 



V3d-J 2 ) 



6 = 2tan~ [ 



1 - 3y 2 + Vl + 3v 2 



3jVT 



P 1 P 
with y:= — = — . 

PL V3Poo 



(3.31) 

(3.32) 
(3.33) 



From them one can plot the renormalized momentum -P/pco = 
y/3yN + 6 and can show that —P p < always holds, i.e., the 
soliton is always stable. One can also obtain the coefficient of 
e 2 in Eq. ( 13.8b as follows: 



d + d\P 



1 



y tanh y 



i£p = 

aP p 4 -y/^p^ y + 2(1 - y2) taml -! y 

V3d-y 2 ) 

7 ■= ^ ■ 



(3.34) 



(3.35) 



The energy-dependence of the reflection coefficient R of lin- 



earized excitations is obtained by solving the Bogoliubov equa- 
tion P. 171 ) numerically, and the results are shown in Fig. [3] We 
can verify that the expression (13.81 ) with ( 13.341 1 is valid for low- 
energy region. From this figure we can also see that the soliton 
with zero velocity is the strongest scatterer. It is intuitively clear 
since the shape of the soliton becomes shallower and wider if 
the velocity of the soliton increases. However, this intuitive un- 
derstanding is not always correct, as the integrable cubic case 
in Subsec. 13.2.21 and the instability-induced anomaly in Sub- 
sec. 13.2.41 illu strate . 



3.2.4. The case with a\ < and > 

Finally, we consider the case with a\ < and aj > 0, which 
is most interesting from the viewpoint of critical phenomena, 
since the soliton can become unstable and the reflection coeffi- 
cient can show the singular and scaling behavior. 

If both a\ and a-i are nonzero, we can set \a\\ = fal = 1 with- 
out loss of generality by the following scale transformation: 



R 



t 

Vk>T 
Id I 



(f)(x, t) — —(p(x, t) 



(3.36) 
(3.37) 



So we performed numerical calculations by setting (01,02) = 
(—1,1). Note that p^ cannot be normalized to be unity if we 
choose £ and rj as the above. Another choice of rj is possible 
to normalize p^ = 1, but in this case either a\ or ai cannot be 
normalized. 

Using the dark soliton solution ( 13.191 ), we numerically solved 
the stationary Bogoliubov equation ( 13.17I ). and constructed the 
solution with the asymptotic form (12.35b . Figure [4] shows the 
reflection coefficient with = 0.45 for various soliton veloc- 
ities. We can observe that the zero-energy phonon transmits 
perfectly, unless the soliton velocity is equal to the critical one. 
When the soliton velocity comes close to the critical value, the 
slope of the reflection coefficient becomes very steep, and at the 
critical state, the perfect transmission eventually vanishes. The 
approximate expression (13.8b is good for sufficiently low exci- 
tation energy. Figure [5] shows the scaling behavior of reflection 
coefficient R. Here, based on Eq. (13.13b . the horizontal axis is 
chosen to be the scaled energy e = e/ ^2P"(p c )(P(p) - P{p c )). 
We can see that if the soliton velocity is close to the critical 
one, numerically calculated points are well fitted to the univer- 
sal curve ( 13.10b . Thus, the theoretical results are well confirmed 
in this example. 



4. Proof - Step 1: Exact Zero-Energy Solutions 




• p = p c = 0.2066 p L 
o p= 1.05 p c = 0.2169 p L 
□ p= 1.1 p c = 0.2273 p L 
t p = 0.25 p L 
a p= 0.3 p L 
3 S ♦ P= 0.4 p L 

%fig . P = 0.5 PL 



"""""""""""""■?'"M-M^ 



Figure 4: (Color online) Energy-dependence of reflection coefficient R of lin- 
earized excitations for various soliton velocities. Here we set (01,02) = ( — 1. 1) 
and poo = 0.45. The critical velocity of the dark soliton is given by p c = 
(0.206597 . . . ) X pl- (See the lowest curve of Fig. [2]) A reflection coefficient of 
zero-energy phonon at the critical state is given by (d/b) 2 =s 0.5718. Parabolic 
curves represent the theoretical approximate expression 13.81 . 



4.1. Parameter derivative 

As seen in the asymptotic form ( 12.21b or in the example 
of the CQNLS system in Subsec. 13.21 the dark soliton solu- 
tion has two parameters, i.e., (/7,pco). So we can consider two 
kinds of parameter derivatives^] d p (f> and <9 Poo 0. We can use 
arbitrary coordinates to "label" the two-dimensional parame- 
ter space (a,/3) = (a(p,pco),/3(p,p<x,)), unless the Jacobian of 
coordinate transformation is singular. Obviously, the final re- 
sult must not depend on the choice of coordinates. In order to 
make the story general, we always use these general parameter 
derivatives, and henceforth, we write the parameter derivative 
simply by the subscript, i.e., (f> a :- d a <p and <f>p := dp(f>. We also 
introduce the following symbol: 



[A, B] a a :- A a Ba - ArB . 



Note that the ratio 



[C,DU 



(4.1) 



(4.2) 



has a coordinate-free meaning, in other words, it is invariant 
under coordinate transformations of parameter space. We often 
construct coordinate-free solutions in such a ratio form. 

An immediate application of parameter derivatives is that 
one can obtain a particular solution of zero-energy Bogoliubov 
equation ( 12.9b . By differentiation of the stationary NLS eq. 
J2.8b . one obtains 



A. 



(4.3) 



In this and the next section, we prove the main result. This 
section is particularly devoted to the construction of exact zero- 
energy solutions. As an important tool, parameter derivatives 
are introduced. 



Here we do not consider the derivative with respect to the parameters in- 
cluded in the definition of the nonlinear term, e.g., a\ and ai in the CQNLS 
equation, because these parameter-derivatives do not yield a solution of the lin- 
earized equation. 
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R 

0.6 r 




Figure 5: (Color online) Scaling behavior of reflection coefficient R. Here we 
set(ai,o 2 ) = (-1,1) and poo = 0.45. e := e/ J2P"(p c )(P(p) - P(p c )) is a 
scaled excitation energy. "Theory" represents the universal form of reflection 
coefficient <3.10t . 



The same expression follows by replacement a — > ft. Thus, 
taking the difference of double parameter derivatives, one can 
obtain the following zero-energy solution: 



lv, <p] a p 



= 0. 



(4.4) 



It must be emphasized that this solution exists even when a 
localized potential barrier is added, in other words, when the 
fundamental equation loses a translational symmetry. What we 
only need is two kinds of parameter derivatives. So, this solu- 
tion is not a symmetry-originated zero-mode. (For a symmetry 



consideration, see Appendix D ) 

Some technical (but crucially important) identities are de- 
rived in Appendix A Equation ( 14.31 ) will be used again in the 
process of energy expansions. 

4.2. Density fluctuation and phase fluctuation 

Here we introduce notations for the linearized density fluc- 
tuations and phase fluctuations, and rewrite the Bogoliubov 
equation with respect to these variables. They are convenient 
for both calculations and physical interpretations. Through the 
symbols (u, v) = (5<p, -8<p*), the density and phase fluctuations 
are expressed as 



8p = 8(<p(j)*) = 5(f>(f>* + 



Therefore, if one defines / and g as 



1 / u v 
2i'- + - 



(4.5) 
(4.6) 



then / has the meaning of the phase fluctuation, and g the den- 
sity fluctuation. The stationary Bogoliubov equation ( 12.91 ) are 
rewritten as follows: 



pf + jj = Jg, (4-8) 



j P 



- ^-g - 2pF'(p)g - 4jf = -liepf. (4.9) 
P 



For the equation with e = 0, we obtain the following integration 
constant: 



pf + 11= const. 



(4.10) 



4.3. Exact zero-energy solutions 

Let us derive all four linearly independent solutions for the 
zero-energy Bogoliubov equation (12.9b . or equivalently, Eqs. 
( 14.81 ) and ( 14.91 ) with e = 0. From global phase symmetry 
and translational symmetry, we can immediately find two zero- 
modes: (u, v) = (i<p, i<p*) and (<p x , —</>*). In addition to them, we 
already have the third solution in Subsec. 14.11 It is well known 
that if we have n — 1 linearly independent solutions for an n-th 
order linear differential equation, the last one can be obtained 
by reduction of order. So, the fourth solution also can be found. 
See | Appendix B| for a detailed calculation. 

Thus, all four zero-energy solutions are given by 



a- 

i 



[p,S] a /3 



[P, j]a{S \IM>PU 
S x- P\ _ (p(Pca-p)/p 

Px ) \ P. 

Pg2 /0' 



P 



r 

Jo ( 



(Pco - pf 



(4.11) 
(4.12) 
(4.13) 
(4.14) 



For reference, we also write down the same solutions by (u, v) 
notation: 



[pj] 



ap 



g2 



[p, 4>]afi 
-[P,$*]c 



2(pco - P) 



-4> 



f 

Jo 



gidx 



(Pco - pf 



" 3 I 

V3J 



(4.15) 
(4.16) 
(4.17) 
(4.18) 



Here g 2 = [p,p] a p/lp, j] a p = i("iV2 - u 2 Vi). 

(J4, gA) is chosen so that the integration constant ( 14.10b van- 
ishes. So, only (fj, gi) contributes to this constant: 



J ~ 2i \(f> + <f>* 

g = U(f>* - v(p 



U = ^ [f+ 2-p 



(4.7) 



Pfi 



Jg2 _ l 
P 

pf! + J ji=0, 

' P 



(i= 1,3, and 4.). 



(4.19) 
(4.20) 
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If we assume that the density p(x) is an even function, the 
parities and asymptotic behaviors of the four solutions are sum- 
marized as follows: 



(A. si) 

(fl,g2) 
(f3,g3) 
(f4,g4) 



(even, odd), bounded, 

(odd, even), linearly divergent, 

(even, odd), exponentially decreasing, 

(odd, even), exponentially increasing. 



(4.21) 
(4.22) 
(4.23) 
(4.24) 



Keeping in mind the above (l4.2H - d4.24t will help to un- 
derstand and imagine the construction of the solution in the 
next section. More detailed forms of asymptotes are given in 



Appendix C 



Deriving all solutions including divergent ones might seem 
to be only of academic interest, but they will become necessary 
for the exact energy expansion, which will be performed in the 
next section. 

We note that four solutions are also obtained without using 



the parameter derivatives as shown in Appendix F However, 
those expressions are almost useless for solving a scattering 
problem. 

5. Proof - Step 2: Low-Energy Expansion 

In this section, with the use of zero-energy solutions derived 
in the preceding section, we construct finite-energy solutions up 
to second order by an exact energy expansion method. Using 
them, we solve the scattering problem of linearized excitations, 
and prove the main result. 

5.7. Exact low-energy expansion 

We construct the finite energy solution by energy expansion: 



I n=0 v 7 



(5.1) 



or equivalently, (/, g) - Y*7=o e " 8^)- The expansion only 
using asymptotic forms was used in Ref. [31], and the com- 
pletely exact expansion was first used in Ref. [32] to discuss 
tunneling behaviors of Bogoliubov phonons in the presence 
of potential walls. In order to prove the singular behavior at 
the critical point, the exact expansion is essentially necessary 

The substitution of (15.11 ) into Bogoliubov equation ( 12.91 ) 
yields 



( (n-l) 
,(«-!) 



(n= 1,2,...), 



(5.2) 



or in (J, g) expression, 



pM'+^Mr, (5-3) 



_ L g inj _ 2 pF'(p)g ( " ) - 4;(/W)' = -2ip/<- 



(5.4) 



It is equivalent to the zero-energy Bogoliubov equation with an 



inhomogeneous term (li 



(n-\) v (n-l) 



'). Since we already know all 



four solutions for the homogeneous equation, we can solve it by 
variation of parameters. Since we want the solution of the form 
( 12.351 ). in the process of energy expansions, we must cancel out 
the exponentially divergent terms, and must take up the power- 
law divergent terms, because e = 1 + i/c, jc - k 2 x 2 /2 + • • • and 
k\, &2 oc e for low-energy. (See eqs. ( 12.171 ) and ( 12.181 ).) 

5.2. First order solutions 

The first order solutions are, in fact, directly calculated with- 
out the use of variation of parameters. When we choose 



(w 1 



(0) v (0) 



; ) = («i,vi) = (i<f>,i<f>*), Eq. ( 14.31 ) gives the particular 



solution for (u (v \ v (1) ), that is, 



and 



~ (finite e) 
M A 

-(finite e) 
V A 



-(finite e) 
B 

-(finite e) 



£♦£1 "!♦<*"> 



(5.5) 



16 



Ms 



+ 0(e 2 ) 



(5.6) 



give two kinds of first order solutions of the finite energy so- 
lution. From Galilean symmetry of the NLS equation, we can 
also obtain the first order solution when we set (w (0) , v (0) ) = 



(4> x , -0*). (See Appendix D I It is given by 

-(finite e) 



(finite e) 



4>, 



We also write down them in (/, g) notation: 



infinite e) 

J A 
-(finite e) 



0/ j_l a \p a ) 



infinite e) 
.(finite e) 



infinite 6) 
^3 

-(finite e) 

«3 



Me \Pp) 



(5.7) 

(5.8) 
(5.9) 
(5.10) 



5. 3. Second order calculation - Identification of bounded solu- 
tions 

Let us calculate the second order. Before beginning, we ex- 
plain the outline of the calculation. It is important to notice 
that only two of four solutions are the bounded solutions, i.e., 
propagating waves for finite positive energy e, as seen in Sub- 
sec. 12.21 This means, in the energy expansion, that the two of 
four solutions are power-law divergent and the remaining two 
must be exponentially divergent. On the other hand, however, 
we have three linearly divergent solutions ( 15.8I )- (15.10I ) in the 
previous subsection. So, one of the three must be a dummy. 
In fact, the solutions calculated from the variation of parame- 
ters shows that all three solutions diverge exponentially; The 
parity of the second-order exponential divergent term becomes 
(/ (2) , g (2) ) = (even, odd), so it cannot be canceled out by adding 



9 



the homogeneous solution (fi, #4). In order to kill out this expo- 
nential divergence, we must choose the special linear combina- 
tions of the three. Through this process, we have two power-law 
divergent solutions, as expected. 

Now, let us construct the solution of Eqs. ( 15.3b and ( 15.41 ) 
with n — 2 by the above-mentioned method. In accordance 
with the method of variation of parameters, let the particular 
solution be (/ (2) ,g (2) ) = 2i=i,2,3,4 c <C/j"» and let c,s have the 
x-dependence. After a little calculation, the equations for c,s 
are summarized as follows: 



If we write 



2ih 



2ip 



h 


S2\ 


l8 il) 


-/1 


gl J 




U 






-h 


83) 





(5.11) 
(5.12) 



Let us assume that is a linearly divergent odd function 
and g (1) is a bounded even function. From the asymptotic 
behaviors (I4.21b - d4.24b and Appendix C we can see that 
c\(fi, gi), C2C/2, 82), and cj,(fi,,g3) are always power-law diver- 
gent functions. C4 is an asymptotically constant function, so the 
product cn(ji, gn) diverges exponentially in general. Exception- 
ally, if C4 decays exponentially, c\(f\, g\) becomes a power-law 
divergent function. Therefore, in order for the particular solu- 
tion to be a power-law divergent function, lim^i-.o Q = 0, or 
equivalently, 

<dx = <-> J (/ 3 g (I) -g 3 / (1) )d* = (5.13) 
must hold. 

Let us construct the solution that satisfies the above condition 
( 15.13b . Consider the following linear combination: 



y^finite e) 
(finite e) I + 2^A I ."(finite e) 



Xfinite e) 

h 



8a 

1 + 2f A S 

£a which satisfies Eq. ( 15.13b is determined as 



S3 



(5.14) 



f 

»y -o 



<->£a = 



[/3Pff - 8i(S a + liaijxx)] dx = 
JV/z„ 



(5.15) 
(5.16) 



The same calculation can be performed with respect to the other 
parameter /?, so (/^ finite e) , g^ finite e) ) and £ B are defined in the 
same manner. 

Thus, we have two non-exponentially-divergent solutions 

/ /-(finite e) (finite j / /-(finite e) (finite eK *.i 

(J a '8a ) and (/g ,g B ). Furthermore, we can 
get coordinate-free solutions by making the linear combinations 
of these two. They are defined as follows: 



(finite e)\ 

(finite e) I 
\Sl / 

^(finite e )v 
^(finite.) J 



1 



1 



2(fA"ft) 



y(finite e) 
^A I ^(finite e) 



/■(finite 6) 

Ja 

(finite 6) 



/■(finite e) 

J A 
(finite e) 

•A 

■(finite e) 



(finite 6) 



(5.17) 

/■(finite e)\ 
I „(finite e) I • 

' (5.18) 



■(finite e) 
(finite 6) 



./; 



°° I An)\ 
«=0 7 



(i = 1 or 3), 



then the zeroth and first order terms are given by 



■(()) 



/ 3 (I) 

8? 



i 



[Poo [5, 5 ]^ + p[JV, S V 

p<x,K/i] aj s + p[N,n] a p \P°o[6,p\ a f} + p[N,p] a/3 1 



83 

ix/1 



S x - P 

Px 



(5.19) 
(5.20) 

')■ 

(5.21) 
(5.22) 



iN[M> flap 



2\0) 2{poo[S,ii\ ap + p[N,p] afi ) \g2) \g 



(5.23) 



and the second order terms are, from Eqs. ( 15.1 lb and ( 15.12b . 
given by 



f, 



(5.24) 



with ; = 1 or 3. 

Every bounded solution for finite positive energy e must be 
constructed as a linear combination o/(/j finlte € \g\ ^ e ') and 

^(finite ^ ^(finite e)^ &Q mufJt ^ solution of me scattering prob . 

lem ( 12.35b . Our remaining work is to calculate their asymptotic 
behavior. 

5.4. Second order calculation - Asymptotics 

The evaluation of asymptotic form of the second order term 
d5.24b is tedious but straightforward. All calculations can 
be carried out by using the expressions in | Appendix C| For 
brevity, we introduce the following symbols for the asymptotic 
forms of f> and L g. dx (i — 1 or 3.): 



jf^ -> i(l n x + lio sgnx), 
'dx -> i(m n x + m i0 sgnx). 



f 

Jo 



(5.25) 
(5.26) 



Explicit expressions for l,js and m,jS are given in Appendix E 
The asymptotic form up to second order is given by 

(finite e)\ 



f. 



(finite e) 
7 

e 

2^2 



Sn\ . Ih\x+ kosgnx 
oj ma 



]t[( p l mn ~ + (Pl" 1 ^ ~ jho)\x\]+Cio 

(pooki - pm n )x + (poJio - pm i0 ) sgn x / 



+ 0(£ 3 ). 

(5.27) 
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Here c,o is a certain constant whose form is not important here. 

Let (uf nite e \ vf nite e) ) be the counterpart of (/f nite e) , gf nite £) ) 
in (w, v) notation. If we write them in plane wave form 



with 



(finite e)\ 
(finite e) j 



and 



(finite e)\ 
(finite e) 



i yfp^ 

2PPL 



i y/p^ 
2PPL 



[Cfwi(jC,±|) + Cjw2(*,±f)] (x^+oo) 

(5.28) 



[0fwi(*±f) + Dfwa(x,±f)] (jC-»±oo), 

(5.29) 



a = 4p L poo, 

b = (N + pN p f + (/> L Ag 2 , 

d= (N + pN p f - ( PL N P ) 2 , 

It gives the result in Subsec. 13.11 

6. Discussions and Concluding Remarks 



(5.39) 
(5.40) 

(5.41) 

(5.42) 
(5.43) 



then the energy dependence of coefficients can be obtained as 
follows: 

Cf = p(p + p h + 2k 2 l n ) ± iepdp - pOho + 0{e 2 ), (5.30) 
Ct = P(~P + Pl- 2k 2 1 ix ) ± ie Ph (p + pOho + 0{e\ (5.31) 

Df = 2p K 2 h, ± iep L + (p - p L )/ 30 1 + 0(e 2 ) 



4p, 



(5.32) 



= -2^/ 3] ± iep L + (p + p L )Z 30 j + 0{e l ). (5.33) 

The striking feature is that the zeroth order of and D 2 van- 
ishes when the soliton velocity becomes the critical value, see 
Eq. (IE.23I ). This is an immediate cause of the singular behavior 
of the reflection coefficient. 

A solution of the scattering problem (12.35b is constructed as 
follows: 



/ (finite e) 
D 2 ^(finite e) 



/ (finite e) 
2 y (finite e) 



(5.34) 



and coefficients t and r are given by 



t = 



C+D- 
C+D--D+C- 



D+C- 

D 2 C 2 



(5.35) 
(5.36) 



Finally, moving to the particular coordinate system (a, ft) = 
(p,Poa), and using the expressions given in Appendix E we get 
the main result 



2f&6 f 
" P 2 Pt 



(C 2 D\ - DjCJ") = aP p - ie(b + biP p ) + 0(e 2 ) 

(5.37) 



In this last section, we give a few discussions and future per- 
spectives. 

6.1. Local density fluctuation at the critical point 

In the system with a potential wall, the emergence of a zero- 
energy local density fluctuation was a key of the destabilization 
of superflow B 3 2l 1 3 311 . In the present case of solitons, an amplifi- 
cation of the zero-energy local density fluctuation also occurs at 
the critical point, but its mathematical structure slightly differs. 
Let us see it in detail. 

Before beginning, one should recall that / and g have the 
meaning of phase and density fluctuations. (See subsection !4.2l 
again.) So, the solution (u, v) cc (<p,<p*) <-» (/, g) cc (1,0) has 
no density fluctuation. If another solution, e.g., a parameter- 
derivative solution, is added, a non-zero density fluctuation 
arises. 

For a superflow state against the potential barrier| 17, 18, "l9h . 
it is knownJUHl] that 



limh-'* 

e->0 \ V 



lim 

E^O \ v 



+ c- 



dtp \-<f>* 



(for non-critical states). (6.1) 
(for a critical state). (6.2) 



Here tp is a Josephson phase difference and c is a certain con- 
stant. Thus, the density fluctuation represents the anomaly of 
the critical point. 

In the case of solitons, however, because of spontaneous 
translational symmetry breaking, the density fluctuating zero- 
mode (fs,gi) (or equivalently, (u^, V3) ) always exists, and this 
mode indeed contributes to the solution of the scattering prob- 
lem: 



lim 



= Ci 



+ C 3 



(for non critical states). (6.3) 



and 



" ; 2 ^2^2 
P 2 Pi 



(C 2 D 2 - D + 2 C 2 ) = -ie(d + dyP p ) + 0(e 2 ) (5.38) 



Here c\ = lim^oi-C^) and C3 = lim e ^o D 2 are constants. (See 
Subsec. [5~4l for more detailed expressions.) Therefore, the local 
density fluctuation always exists regardless of whether the soli- 
ton is stable or unstable. However, when the soliton velocity 
comes closer to the critical one, the ratio C3/C1 becomes larger, 
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and it becomes infinite at the critical point. That is to say, at the 
critical velocity state, 



lim 

mo \v 



(for a critical state) 



(6.4) 



holds. Thus, we can say that the amplification of the local den- 
sity fluctuation also plays a key role in the case of the destabi- 
lization of solitons. 

6.2. Conclusions and Future perspectives 

In this paper, we have solved the scattering problem of 
linearized excitations (Bogoliubov phonons) against a dark 
soliton in a generalized NLS system. We have exactly shown 
that the perfect transmission of a zero-energy phonon vanishes 
when the soliton velocity reaches the critical value, and near 
the critical velocity state, the reflection coefficient obeys a 
saddle-node type universal scaling law. Our result has a 
fundamental importance because it provides an exact example 
of saddle-node scaling in infinite dimensional time-reversible 
Hamiltonian systems. Through the proof, we have also 
obtained the exact zero-energy solutions and their finite energy 
generalizations. In the derivation of them, the use of two 
kinds of parameter derivatives has played an important role. 
This method will also be useful to elucidate the low-energy 
physics of other systems, such as higher dimensional systems 
or multi-component systems. 

Even though we have shown the example of scaling laws, the 
derivation of a normal form of saddle-node bifurcation remains 
unsolved. The similar problem for the supercurrent-flowing 
system in the presence of an obstacle lll8l also exists. 
Compared to the problem of solitons treated in this paper, the 
system with an obstacle seems to be a little more difficult, 
because all four zero-energy solutions for linearized equation 
are not yet obtained except for the critical velocity state[32]. 
These issues are left as future works. 

In this paper we have treated the soliton-phonon scattering 
problem. One interesting generalization is the multi-soliton 
scattering process. It was shown that the existence of the two 
solitons makes the reflection coefficient of phonons non-trivial 
even for the integrable cubic case ll34n . The soliton collision 
problem in the non-integrable generalized KdV equation was 
recently investigated 11351 13611 . The similar problem of dark 
solitons in non-integrable NLS systems is also important to 
understand the generic solitary characters of dark solitons. 

An emergence of singularity of scattering properties at a 
critical state which separates the stable and unstable branches 
is expected to be a universal character in more general systems, 
because the emergent or amplified zero-modes can affect 
the transmission properties of low-energy modes. Since the 
scattering problem of linearized excitations is easier and more 
analytically tractable than the existence-proof of an unstable 
mode or construction of a Lyapunov function, it will be useful 
to "conjecture" the criterion for stability of solitons, even 
though it is not a direct proof of the stability itself. For exam- 
ple, to the best knowledge of the present author, the stability 



criterion of dark solitons in multi-component NLS systems 
is an untouched problem. In ultracold atoms, the binary 
mixture of Rb atoms is realized l37i l38l] and it is known that 
the corresponding coupled NLS equation becomes integrable 
only for the Manakov case. Also, the Bose condensates with 
spin degree of freedom are created by optical trap[39], and 
Wadati group members have studied the spin-1 solitons at 
the integrable point] 40, 41, 42]. When the coupling constant 
deviates from the integrable one, it is no longer ensured that 
a dark soliton is always stable. So, it may be an interesting 
problem to study the stability of solitons through the study of 
scattering problems of linearized excitations. 
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Appendix A. Identities for parameter derivatives 

Here we derive a few identities for parameter derivatives. 
The NLS equation ( 12.81 1 expressed in terms of the density p(x) 
is given by 



2 -2 
Px . J 



(A.l) 



From the parameter derivative of Eq. ( 12.281 1 and the above equa- 
tion, one obtains 

PaPxx ~ PaxPx = 2(p ra - p)(2pj a - p a p). (A.2) 

Here recall that j = p^p and p — p 2 + F(pco). The same expres- 
sion follows by replacement a — > /3. From them, an important 
identity 



8382 - 8283 = -4p(Pco - p) 



(A.3) 



follows, where g2 and gi are defined in Subsec. 1431 

Dividing both sides of dA.21 ) by p, and integrating them from 
— oo to +oo, we obtain 



2j a 6 + p a N ■■ 



d_ r 00 (p,) 2 dx 

da J_ M 2p 



(A.4) 



Since the same expression for /3 holds, we obtain the second 
important identity 



(2j a 5 + p a N)n = (2j p 5 + p p N) a 
<-> 2[j,5] a p + [p,N] a/ 3 = 0. 



(A.5) 



Appendix B. Derivation of the fourth zero-energy solution 

Let us derive the expression for the fourth zero-energy so- 
lution (f4,g4). Eliminating / from Eqs. (14.81 ) and ( 14.91 1, one 
obtains the third order differential equation for g: 



g"'+Ag' + B = 0. 



(B.l) 
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Here A and B are some functions, but their forms are not im- 
portant here. Knowing two solutions g2 and #3, the last one is 
obtained by reduction of order: 



84 := 82 — — - g3 

Jo w Jo 



with 



vf = 8382 ~ 8283- 



(B.2) 



(B.3) 



Since w is given by (IA.3I ). the above expression is simplified as 



16p 2 g 4 



1 



1 



Poo - P poo - PO 



82-83 



82<lx 

(Poo-P) 2 



(B.4) 



Here po := p(0). The solution which cancels out the constant 
J4.10I ) is given by the following linear combination with g2' 



84 :-- 



Po 



•Po 



-82 - 16p Poo§4, 



(B.5) 



and /4 is calculated from f. = -j84/p 2 



Appendix C. Asymptotics of zero-energy solutions 

Let us write the asymptotic form of p(x) as 

p(i)=p to -pi! ) e- 2 ' w + ..., (C.l) 
then a half of the healing length k is given by 

(C2) 



From d2.21l l or 1 12.271) . an asymptotic form of the phase S (x) is 
given by 



S(x) -> px+ (sgax)-. 
Using the invariant nature of the ratio 1 14.21 ). one can show 



(C.3) 

(C4) 
(C5) 



With the use of the above relations, the asymptotic forms of 
zero energy solutions are evaluated as follows: 



IP, P]af) 


Im, p] PP „ 


Pi 


[p, jlafi 


[p,j]pp x 


Poo a: 2 


lP,Pco] a p 


lP,Poc]pp» 


p 


[p, j] a p 




K 2 



h 



82 



Pi / 1 foflq 

— -x + (sgnx) 



Poof 

_P_ 

K 2 ' 



2[p,j] 



83 -> (sgn x)2Kp y Je 
f 4 -» -(sgnx)^-e 2 * 



g4 



2p^V2 



(C.6) 
(C7) 

(C8) 
(C.9) 
(CIO) 

(C.ll) 



Appendix D. Symmetry consideration on zero-modes 

In this appendix we consider how a symmetry plays a role 
in finding a solution of a linearized equation. Particularly, we 
emphasize that the information we can obtain from a Galilean 
symmetry is not about the zero energy solution but about the 
first order correction of a finite energy solution. 

Let <p(x, t) be a solution of time-dependent NLS equation 
J2.ll ) and ej>(x, t, a) be a family of solutions with continuous pa- 
rameter a such that <f>(x, t, 0) = (p(x, t). Here we assume that a 
is free from any system parameter. (The parameter derivative 
introduced in Subsec. 14.11 is a more general concept, because 
it can depend on the system parameters which appear, e.g., in 
a boundary condition J2.21I ).) Differentiation of NLS equation 
J2.ll ) with respect to a immediately yields 



(id, - £) 



= 



(D.l) 



with a definition <f> a := [d a <j>] a= {). Thus, <p a is a solution of 
a time-dependent linearized equation (12.5b in the presence of 
the condensate wavefunction <j)(x, t). Particularly, if one sets 
(j>(x, t, a) = (p(x, t)e m , (f>(x, t, a) = <p(x + a, t), and Eq. J2.22I ). 
which represent a global phase symmetry, a translational sym- 
metry, and a Galilean symmetry, respectively, one obtains the 
following particular solutions: 



(id, - X) 
(id, - £) 



0, 



(D.2) 
(D.3) 
(D.4) 



It is a result for a time-dependent equation. In order to inter- 
pret these results to a stationary problem, let us set <p(x, t) — 
4>(x)e^' . From the assumption stated above, p does not depend 
on a. (On the other hand, p can depend on the parameter in 
Subsec. 14.11 ) The above equations are rewritten as 



4£W 



4>x 



0, 



(D.5) 
(D.6) 
(D.7) 



Here X.ft defined by Eqs. ( 12.9b and J2.10I ) is a differential oper- 
ator of the stationary Bogoliubov equation. JD.5I ) and JD.6b im- 
mediately give zero-energy solutions, whereas we need a con- 
sideration on Eq. JD.7I ). Since this equation must hold for any 
time t, both coefficients of r and f 1 must vanish. The equation 
for the t l -coefficient reproduces JD.6I ). From the f°-coefficient, 
one obtains 



I r I 



(D.8) 
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It represents Eq. ( EHb with n = 1 and (m (0) , v (0) ) = {<p x , -0*). 
Therefore, it gives the first order solution of ( 15.71 ). Thus, a 
Galilean symmetry gives us a piece of information about first 
order solutions, and gives no new information about zero- 
energy solutions. A linearly divergent zero-energy solution, 
which exists even when the system does not have a transla- 
tional or a Galilean symmetry, can be obtained by two kinds 
of parameter-derivatives stated in Subsec. 14.11 

We further note that Eq. (15.7b itself does not describe a phys- 
ically meaningful non-divergent solution. As we show in Sub- 
sec. 15.31 we must construct a linear combination of solutions so 
that the second order term is free from exponential divergence. 
A correct first order solution with non-divergent character is 
given by ( 15.231 ). 



Appendix E. Formulae for calculation of asymptotes 

Let us derive some formulae for ZyS and myS defined by Eqs. 
( 15.251 ) and ( 15.261 ). Using the relations 

[X, j] a p = Poo[X, p] a p + p[X, (E.l) 

[X,n]ap = 2p[X,p] a0 + F'(poa)[X, (E.2) 

one can show 

Poc[X,p] a/} - 2p[X, j] a0 = 2/< 2 [X,p 00 ] a p. (E.3) 

Here remember that p^ and F'(poo) are related to each other by 
( 12.141 ). and k 2 — p 2 ^- p 2 . Using it and dA.51 ). one obtains 

P<x,[8,p] a p + p[N,p] a p - 2k 2 [6, (E.4) 
From (15.21b . (15.23b , and dE.41 ). we can obtain 

p<x>[8,p]afs + p[N,p] afi 



mn 



/io - 



mio = 



2K 2 [6, pcvja/j 
Pco[5,Poo]a/3 + p[N,Poo]a/} 
2K 2 [6,p m ] a/3 

p[N, S\ ap 



4K 2 [6,Poc]al3' 

Poo[6,N] afi 



AK 2 [8,Pao] 



1 N[p,p] ap 

2 Ak 1 [6, poc\ a p 

N[p,Poc]ap 
m 31 = ~ a ~ i /"Mil, 



'30 - _ 



4-K 2 [5,Pco] a p 

N[p, 6] a p 

o 2v x 1 P l W, 



N[M,N] a p 

8K Z [6,Poa]aP 



(E.5) 
(E.6) 
(E.7) 
(E.8) 
(E.9) 
(E.10) 
(E.ll) 
(E.12) 



Further, using them and (IE. II ). (1E.2I ). and ( lE.3b . 

2pl n +F'(p 00 )m n = h (E.l 3) 

Poo/11 +pmn = — , (E.14) 

2 |APooJa0 

1 [N,6U 

2pl m + F'(p x )m w = -- 1 J / , (E.15) 

Poo/10 + pmio = 0, (E.16) 

2p/ 31 + F'(poo)m 3 i = 0, (E.17) 

Poo/31 + pm 3 i = --— — , (E.18) 

2 [O,pcolap 
1 [pN, 6] a e 

2pl 30 + F'(pco)m 30 = - rc : n , (E.19) 



Poo/30 + f"w 30 = — . 



2 [<5,Poo]a^ 



(E.20) 



We can eliminate mys by Eqs. ( IE. 131 ). ( IE. 16b , ( IE. 17b . and 
JE.20b . and the expressions (15.30b - d5.33b are derived in that 
way. 

Now we consider the particular coordinate {a,/3) = (p,p M ). 
Zys are then given by 



hi = - 
ho 



31 



jSp + pjN p 
2p 00 K 2 6 p 
p[N,6] pp „ 
Ak 2 5 p ' 

2p 00 K 2 5 p ' 



N\ti,8\pp„ 



From P = -pN - p^d and (lA.5b . we obtain 
N + pN p + P p 



Poo 

pN - k 2 N p + pPp + jN p , 



(E.21) 
(E.22) 
(E.23) 
(E.24) 

(E.25) 
(E.26) 



Derivatives of 6 are eliminated by using them. With the use of 
Eqs. dB3fl -( lK26l ). we can obtain the final result £37M£43]|. 

Appendix F. Zero-energy solutions without using parame- 
ter derivatives 

We can write down the general zero-energy solution without 
using parameter derivatives: 



:)= z 



1=1,2,3,4 
M[ = l4 
U 2 = <f). 



(F.l) 

(F.2) 
(F.3) 



/dx , C dx A f dx 

— +4i/ 2 — j-4#, -3, (F.4) 
p J pp; J pj 

/dx P pdx 

-2 -0x rV- (F.5) 
p; J p; 
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These expressions, however, have a fatal flaw; Namely, they 
have artificial singularities at the origin (or more precisely, the 
points where p'(x) = 0.) It means that these expressions do not 
give global solutions which continuously connect the solutions 
from x = -oo to x — +00, instead, only give local solutions 
which satisfy the differential equation at each point. For this 
reason they are not so useful in scattering problems, which are 
equivalent to the determination of global behavior of given so- 
lutions. On the other hand, the parameter-derivative solutions 
given in Subsec. 14.31 have no singularities if the density p(x) 
does not cross the value p^. 

Appendix G. Stationary solutions of CQNLS equation 

In the CQNLS system, the equation for momentum conser- 
vation (12.24b becomes 



(A) 2 



= -/ + jmP ~ PP 2 + flip 3 + a 2 p 4 



(G.l) 



Here we want general stationary solutions, so we do not con- 
centrate on the solution with the asymptotic form (I2.21I ). and 
therefore p, j, and j m need not be given by ( 12.1 It . ( 12.251 ). and 
( 12.261 ). respectively. If the right hand side of (IG.ll) is factored 
as fl2 IljLiG _ Pi), the solution of this differential equation is 
given by the following cross-ratio form: 



2, 1 \ (P(x)-pi)(p3 -pi) 

sn (Kx\m) = — — 

(p(x) - pi)(P3 - Pi) 



(G.2) 



with 



k = y/a 2 (p4 - pi)(p3 - Pi), (G.3) 
(P4-pi)(p3 -pi) 



(P4 - pi)(P3 ~ Pi)' 

or equivalently, 

P2O03 -Pi)-pi(p3 -P2)sn 2 (Kx\m) 



p(x) 



(P3 - Pi) - (P3 - pi) sn 2 (Kx\m) 



(G.4) 



(G.5) 



Here we use Mathematica's definition for Jacobi elliptic func- 
tions0 

For the case of the dark soliton solution, p3 = p\ — p^ holds, 
so m = 1 follows and the sn function reduces to the tanh func- 
tion, pi and po '■= P2 are given by the roots of 



a2p + (2a 2 Poo +a\)p- p = 0. 



(G.6) 



po (13.21b is the root of Eq. (IG.6I ) with a plus sign. p\ can be 
eliminated in two ways; From Eq. dG.3b or from the fact that po 
and pi solve Eq. (IG.6I ). 



Poo - pi = 



A2(Poo - Po) 



P0P1 = ■ 

a 2 



(G.7) 



The expression JG.5b is then rewritten as 

K 2 (pco — po) sech 2 kx 



P(X) - Poo - 



k 2 — fl2(p<>o — po) 2 tanh 2 kx 



(G.8) 



The following expressions are also useful for calculation of the 
phase shift 5 and the particle number of soliton N: 

c j , (pof)[p(Pco-po)tanh*-x]' 
5 X = _ =p+ _______ _ ^ — - | (G.9) 



p-Po. 



p (p K ) 2 + [p(poo - po) tanh kx] 2 ' 

1 Kiyfaiipac, - po) tanh kx]' 
y/a2~ k 2 - [ -^(poo - po) tanh kx] 2 



(G.10) 



where Eq. ( lG.9b is obtained by eliminating both a 2 and p\ with 
the use of dG.71). Here we can use the formulae 



ay 



^=(tan-^)', (Gil) 

°y' _L„u-i^ 



a 2 + y 



a 2 - y 2 



(G.l 2) 



Thus the content in Subsec. 13721 is reproduced 



2 The Wolfram Functions Site, http://functions.wolfram.com/ 
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